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Abstract 

We deal with the problem of determining a time varying inclusion 
within a thermal conductor. In particular we study the continuous de- 
pendance of the inclusion from the Dirichlet-to-Neumann map. Under a 
priori regularity assumptions on the unknown defect we establish loga- 
rithmic stability estimates. 
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1 Introduction 

In this paper we study the stability issue for the inverse problem of recovery 
the discontinuous conductivity coefficient of a parabolic equation from infinitely 
many boundary mesurements. 

First let us give a coarse formulation of the problem which we are going to 
study. Let T be a given positive number. Let O be a bounded domain of W 1 , 
n > 2, with a sufficiently smooth boundary and let Q be a domain contained in 
Q x (0, T). Assume that for every r £ (0, T) the intersection D(t) of Q with the 
hyperplane t = r is a nonempty set and f2 \ D{t) is connected and denote by k, 
k =/= 1 a positive constant. Let u be the weak solution to the following parabolic 
initial-boundary value problem 

d t u - div((l + (k - 1)xq)Vw) = infix (0, T), 
u(-,0) = inH, 
u = g on 90x (0,T), 

where g is a prescribed function on dQ x (0,T). The inverse problem we are 
addressing to is to determine the region Q when infinitely many boundary mea- 
surements < g, tt > are available. The problem formulated above arises 

L OV\dQx{0,T)) 

in nondestructive testing evaluation ( |Ca-Moj . |Pa-La-Alj ). 
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A uniqueness result for the problem introduced above has been proven in 
1997 by Elayyan and Isakov [EFTs] . The main tools on which the approach of 
[El-Is] is based are the Runge approximation property and the use of solutions 
with Green's function type singularity. For the nonconstructive character of 
the Runge property, such an argument does not seem suitable for our purpose 
of proving an accurate stability estimate of Q under some a priori informa- 
tion. Thus, along the line of previous elliptic and parabolic inverse problems 
[Al-DC] . [A1- Vej . |Ve| , we abandon such an approach and we choose to use ar- 
guments based on quantitative estimates of unique continuation Al-Be- Ro-Ve] , 
|DC-Ro-Ve] . |Mo-Ro03l IMo-Ro04] . We also make use of singular solutions of 
Green's type, but more quantitative information are necessary in order to obtain 
stability estimates. In particular we need an accurate study of the asymptotic 
behaviour when the singularity gets close to the interface dQ. 

In the present paper we prove that, under mild a priori assumption on the 
topology and the regularity of Q, such an inclusion depends continuously on the 
boundary data with a rate of continuity of logarithmic type (see Theorem l2.7l for 
a the precise statement of the result). In the context of elliptic inverse problems, 
it has been shown that logarithmic stability estimates are optimal ( |DC-Ro| ). For 
parabolic inverse problems with unknown boundaries (and the whole Dirichlet- 
to-Neumann map) examples showing that the continuous dependance can be 
at most of logarithmic type, have been obtained in [DC-Ro-Ve] and [Vej . Their 
proofs work in our situation as well. Namely in such papers the limit cases 
k = +00 and k — are considered. Everything remains basically the same in 
the intermediate situation with < fc < 00, k =/= I. 

A crucial tool to obtain the logarithmic stability estimate is connected with a 
precise evaluation of smallness propagation based on the two-sphere one-cylinder 
inequality for solution to parabolic equations (Es-Fe- Ve] . |Ve] (Theorem 3.10, in 
the present paper). Indeed, roughly speaking, such an inequality allows us to 
approach the boundary of the inclusion in any slice of time. 

Finally we wish to mention here papers of Daido, Kang and Nakamura 
[Da-Ka-Na] and, more recently, Isakov, Kim and Nakamura [Is-Ki-Na] which 
are strictly related to the present one. In [Is-Ki-Na] . the authors consider a 
similar inverse parabolic problem of detecting an inclusion, that does not de- 
pend on the time, by mean of infinitely boundary measurements and provide a 
reconstruction procedure to identify it. 

In this paper we have decided not to deal with the case n = 1. Let us just 
observe that such a case is easier and it can be treated essentially with similar 
assumptions regarding the topology of the set f2\ D(t) (see Remark 12.61 below) . 

The plan of the paper is the following. In Section 2 we state our main result. 
We first give the notations and definitions we need throughout the paper (Sub- 
section 2.1) and then in Subsection 2.2 we state the hypothesis and the stability 
theorem (Theorem 12. 7[) . In Section 3 we provide a proof of Theorem 12.71 We 
derive first some identity which will be the starting point of our proof. Then we 
give some auxiliary result concerning the Hausdorff distance (Proposition 13. 2]) . 
13.31 and 13. 5[) . fundamental solutions (Proposition 13.6] ) and unique continuation 
properties (Theorem 13. 10p . Afterward, using the assumptions on the regularity 
of the inclusion, we derive some further property related to the distance of two 
inclusions and state Proposition l3 . 91 which provides lower bounds for the solution 
of the problem. Finally we prove Theorem 12.71 Proof of auxiliary propositions 
are given in Section 4- Proposition ^ . 31 is proven in Subsection 4.1. In Subsection 
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4-2 we prove Proposition 13.61 and we also give an asymptotic estimate for the 
fundamental solution (Theorem I4.3[) which will be used in the next Subsection 
4-3 for the proof of Proposition 13. 91 

2 The Main Result 

2.1 Notations and Definitions 

For every x E R™, with n > 2, x — (x\, . . . , x n ), we shall set x — (x 1 , x n ), where 
x' = . . . , x n -i) E R n_1 . We shall use X = (x, t) to denote a point in R n+1 , 
where and t e R. For every x € R™ and X = (x, t) E R n+1 , we shall set 

/ n \ 1/2 

M ■ l^l = (N 2 + N) 1/2 - 

Let r be a positive number. For xo E R n wc shall denote B r (xo) — {x E R n : 
\x - Xq\ < r} and B' r (x' ) = {x' E M" _1 : \x' - x' \ < r}. We generally set 
B r = B r (0) and B' r = B' r (0). We denote by B+(0) = {x E B r : x n > 0} and 
B-(0) = {x E B r : x n < 0}. For a point X = (x ,t ) E W l+1 we shall set 
Q r (Xo) = B r (x ) x (t -r 2 ,t ). 

Given a subset A of R", we shall denote by 

[A] e = {xe R" : di$t{x,A) < e} , 

(A) e = {x E A : distfz, dA) > e} , 

[dA] e = {xE R" : dist(x, OA) < e} . 

Let / be an interval of R and let {D(t)} t £i be a family of subsets D(t) of 
R n , we shall denote 

D(I) = \jD(t)x{t} ) 
tei 

and Q = D(R). 

Given a sufficiently smooth function u of variable (x,t) E R ra+1 , we shall 
denote by d l u = df^u = gf.g , i, j = 1, • . . ,n and d t u = ||. For a multi- 
index /3 = (J3i, . . . ,/?„), Pi E N, i — 1, . . . ,n and k E N, we shall denote, as 
usual, dgd^u = gJ^gZ &k , where \j3\ = Yn=i A- Also we shall write V = V x , 
div = div^;. For a matrix A we shall denote by A* the transposed matrix of A. 
We denote by R" = {x E R" : x n > 0}. 

We shall use letters C, Cq,Ci,... to denote constants. The value of the 
constant can change from line to line, but we shall specify the dependancc 
everywhere they appear. Sometimes we have dropped the dependance on n 
which is fixed (n > 2). 

Functional Spaces 

Let D be a subset of R n+1 , / a function defined on D with values in R or R" 
and a E (0, 1]. We shall set 

[fU ° = SUP { (ji f -y\K\t { - 1 (M) ' (y ' s) G D ' ^ (y ' 5) } ' 
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If a e (0, 2] we shall set 



Let fc be a positive integer number, / a sufficiently smooth function and a € 
(0, 1]. We shall denote 

[f]k +a ;D= £ [^/ks, </> fc+Q ;D= £ 'V'V/ ! ...:/>• 

|/3|+2j=fc |/3|+2j=fc-l 

The following Sobolev spaces will be used (we refer to |Li-Ma| for further 
details). We denote by fl a bounded domain in R™. The space H = H q ,3 ' 4 (9f2x 
(0,T)), its dual H' = Hi = H- 3 / 2 '- 3 / 4 {dn x (0,T)), and if = H { / 2 ^/ 4 {dn x 
(0, T)). We consider now the interpolation spaces between H and Hi. For any #, 
< 6 < 1, we define ifg as [iJ , ffi]e, where the latter denotes the interpolation 
at level 9 between the two spaces Hq and H\. The norm in Hg will be denoted 
by || • ||g. First, we notice that for any 0, < 6 < 1, there exists a constant Ce, 
which depends on 9 only, such that the following interpolation inequality holds 
for any ip S Hq 

(2.i) ii^ii^^iHiniv.ii?. 

We also make use of the following notation 

w(n x (o,r)) = {v :v& L 2 ((o,r), ^ e ^((o.r),^ -1 ^))} • 



Boundary Regularity 

Let us give the following definitions 

Definition 2.1 Let be a domain in R™. Given a, a £ (0, 1], we shall say that 
dQ is of class C 1,a with constants po,E > if for any P £ dft, there exists a 
rigid transformation of R n under which we have P = and 

n n B po (0) ={j£ B po (0) : z„ > p(z')}, 

where <p is a C 1:Q function on B' po (0) which satisfies the following conditions 
<f(0) = |V x 'v(0)| = and \\(p\\ C i,*(B> po (o)) < Ep Q . 

Remark 2.2 We have chosen to normalize all norms in such a way that their 
terms are dimensional homogeneous and coincide with the standard definition 
when po = 1. For instance, for any <p E C 1,a (B' po (0)) we set 

IMId.<»(B£ o (0)) = Po|MU°°(B£ (0)) +/ , o||V a! 'V'|U<»(B^ (0)) + Po + ° I V ^' <P\<*;B' po (0) 

Similarly we shall set 

IN™ = po (n+1)/2 (j^dx^'\ 

where <iX = deceit. 
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Definition 2.3 Let Q be a domain in R" +1 . We shall say that Q (or equiva- 
lently dQ) is of class K. with constants po, E if for all Pq E dQ there exists a 
rigid transformation of space coordinates under which we have Po = (0, 0) such 
that 

QH (B po (0) x (-pIpD) = {X E B po (0) x {-pIpD : x n > tptf,t)}, 

where tp is endowed with second derivatives with respect to xi, i = 1, ••• ,n, 
with the t-derivative and with second derivatives with respect to Xi and t and it 
satisfies the following conditions <^(0, 0) = |V X '<^(0, 0)| = and 

Po\\V x '^\\ l ^ { b' po x{-pI p I)) + pU d IM\l^{b' po x{-pIpI)) 

+ Poll%>||L~(i^ o x(-p=,p=)) +Po\\Vx'd t <p\\ LO o {B , pgX{ _ p 2 !p 2 i)) < E Po . 

Definition 2.4 (relative graphs) Let j3 E (0, 1]. We shall say that two bound- 
ed domains fii and fl2 in R n of class C 1 '* with constants Ro,E are relative 
graphs if for any P E dfli there exists a rigid transformation of coordinates 
under which we have P = and there exist <pp,i,(pp.2 £ C 1 ^ (B' ro (0)) ; where 

< 1 depends on E and (3 only, satisfying the following conditions 

Rq 

i) t p Ptl (0) = , |^p >2 (0)|<^ } 

m) \\ l PP,i\\c^e(B' ro (o)) - ER o> i = 1 , 2 > 
in) Qi n B ro (0) = {x E B ro (0) : x n > <p Pyl (x')}, i = l,2. 
We shall denote 

( 2 - 2 ) 7(fii,Q 2 )= sup \\<fip,i- <pp,2\\ LOO ( B , (0) \ ■ 



The Dirichlet to Neumann map 

For any g E H , let u E W{Vl x (0, T)) be the weak solution to the initial- 
boundary value problem 

(2.3a) dtu - div((l + (k - l)xo)Vit) = infix(Q,T), 

(2.3b) u(x, 0) = for x E fi, 

(2.3c) u(x, t) = g(x, t) for (x, t) G dSl x (0, T), 

where \Q is the characteristic function of the set Q. 
Then, for any g E H , we set 

du 

(2.4) A Q g = g^\du X [ ,T), u solution to (JOJ. 

We have that there exists a unique solution u E W(Q x (0, T)) to problem (|2.3p 
[Ev| . In addition, by standard regularity theorems [Li] , [La-So-Ur] and by trace 
theorem jLi-Ma[ Chapter 4, Theorem 2.1], we conclude that AQg belongs to the 
space H and that the operator Aq : H — > H is bounded. We can also consider 
Aq as a linear and bounded operator between H and H' = Hi, by setting 
(2.5) 

{^q9,4>)h',h = (-^-\dnx(o,T),4>)H',H = / -£-<t>, for any g, E H, 

av JdQx(0,T) av 
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where u solves (|2.3[) and (•, -)ff',H is the duality pairing between H' and H . 

Let us remark that the operator Aq is usually referred to as the Dirichlet- 
to-Neumann map associated to the equation (|2.3a[) . 



2.2 Assumptions and Statement of the Main Result 
Assumptions on the domain 

Let po, M, E be given positive numbers. We assume that Q is a bounded domain 
in R ra satisfying 

(2.6a) |fi| < Mp%, 

where |f2| denotes the Lebesgue measure of f2. We also assume that 
(2.6b) <9f2 is of class C ' with constants po, E. 

A priori information on the inclusion 

Denoting by Q = UteR x W (Q = D((—oo, +oo))), we assume the following 
conditions 

(2.7a) dQ is of class K, with constans po, E, 



(2.7b) dist(D(*),9fi) > p , D(t) C ft, Vte[0,T], 



(2.7c) Q,\D(t) is connected Vi € [0, T]. 

Remark 2.5 Let t be any number in [0,T]. Observe that (|2.7a[) automatically 
implies a lower bound on the diameter of every connected component of D(t) 
and Q\D(t). In addition, combining (|2.6a|) with (|2.6bj) . we have an upper bound 
of the diameter of O and thus of D{t). Note also that (|2.7a|) and (|2.7bjl implicitly 
comprise an a priori upper bound on the number of connected components of 
D(t). 

Remark 2.6 For n = 1, it is possible to obtain Theorem 1 2 . 71 replacing assump- 
tion by considering Q = (Q,L) and D(t) = {x e K : s x (t) < x < s 2 (t)}, 
where Sj, i = 1, 2, are C 1 functions such that for all t € (0, T) 

L-s 2 (t)>p Q , si(t)>p a , s 2 (t) - si(t) > po 

and 

\\ s i\\L*°((0,T)) + Poll s 'ilU~((0,T)) < Ep , i = 1,2, 

here 3/c>o < L. 

Theorem 2.7 Let fi C K™ satisfying (|23|) . Let fc > 0, fc 7^ 1 &e c/wen. Let 
{L*i(t)}tgR, {-D2(t)}tGR be two families of domains satisfying (|2.7|) . issame 
t/iat /or e > 0, 



(2-8) IIAQi- Aq 2 \\c(h,w) <e, 

where Qi = Z?i((— 00, +00)), i = 1,2. TTien 

(2.9) d n (D^),D 2 j£j)<u Jt (e), te(0,T], 
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where Lu t (s) is such that 

(2.10) u t {s) < Cpollogs]-", 0<s<l, 

with C = C(t), depending on t,M,E,k only, and < r\ < 1 depending on 
M, E, k only. In addition we have that C(t) tends to +oo as t tends to 0. 

Here d-n denotes the Hausdorff distance. 



3 Proof of Theorem [277 



For the sake of brevity we name a } ■ = 1 + (fc — 1)xq 3 = 1 ; 2. We fix g E H. 
We shall denote by Uj, j — 1,2 the solution of (|2.3[) when Q = Qj. For ^ S 
ff 1 ' 1 ^ x (0,T)) such that 

(3.1) V(- I 7 1 ) = inO, 

using the weak formulation of (|2.3[) we have 



aj-—^-ipdS+ I Uj(x,0)ijj(x,0)dx 



aox(o,T) ^ Jf! 



(ttj Vztj • V?/> — Ujdtip) dxdt = for j = 1, 2. 

Qx(0,T) 



Subtracting the two equations we obtain 

(3.2) / (oiV(mi — u 2 ) • Vip — (ui — U2)d t tp) dxdt 

i!2x(0,T) 



(oi - a 2 )Vw 2 • =< (A Ql - A Q2 )g,tp >h',H, 

Qx(0,T) 

(we notice here that in these identities it is possible to have Uj (•,()) 7^ for 
i = 1,2). Taking ip such that it satisfies ()3.ip and 

(3.3) d t tp + div(aiV^) = infix(Q,T), 

by (|3.2p we have (recalling that on dft x (0, T) u% = u 2 = g) 

/ (ai-a 2 )Vu2-^=<(AQ 1 -AQ 2 )9,ip>H',H, V 5 e H 
Jnx(o,r) 

or, equivalently, 

(3-4) / / (xQi - XQ 2 )^ U 2 ■ V^dxdt = —— < (A Ql - A Q2 )u 2 ,tp >h>,h ■ 

J «/ Si 

Let us denote by T 2 (x, t; y, s) and rf (x, t; y, s) the fundamental solutions of 
the operator <9 t — div(a 2 V) and dt + div(aiV) respectively (rj (x, t; y, s) — if 
t > s and r 2 (x, t; y, s) = if t < s), that is 

[-r 2 (ar, t: y, s)d t <t>{x, t) + a 2 V x Y 2 {x, t: y, s) ■ V x (f>(x, *)] dxdt = cj>(y, s), 
•i+i 

[T{(x,t;y, s)d t <j)(x,t) + aiV x ri(x,t; y, s) ■ V x (j){x,t)} dxdt = <j>(y,s), 
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for every <j> € Cq(M™ +1 ), that is using the 5 Dirac symbol, we have respectively 
d t T 2 (x, t; y, s) - div(a 2 V x r 2 (a;, i; y, s)) = 5(x -y,t- s) 

and 

-d t Tl(x,t;y,s) - div(aiV a: r*(£c,f;y, a)) = 5(x - y,t - s). 

Choosing in (|3.4[) u 2 (x,t) = T 2 (x,t;y,s) and ip(x,t) = T*(x,t;£, r), with (y,s) 
and (£, r) ^ f2 x (0, T), < s < t < T, we obtain 



( 3 - 5 ) / / (XQi - XQ 2 ) v ^ r 2(a;,*;y,s) • V*r*(a;,t;£,T)da;£ft 
Jo Jo 

= t— ^ < ( A Qi ~ A Q2) r 2(-,-;y,s),ri(-,-;C I ' r ) 



For t £ [0, T] we shall define (?(t) as the connected component of fl \ {D\{t) U 
D 2 (t)) that contains 90, Q{t) = (R™ \ O) U C/(i) and 0((O, T)) := Ute(0,T) x 
{*}. For (y, s), (£, r) G £((0, T)) with < s < r < T, we set 

: = / y x T 2 {x,t;y,s)-V x r* 1 (x,t;^,T)dxdt, 



/ VaA^ijy, s) ■ V x Tl(x,t;£,r)dxdt 
JQ2 



U{y,s;£,r) := Si(i/,s;{,t) - S 2 (y,s;£,T). 



Remark 3.1 Let us observe here that for t < s, Si and S 2 can be defined 
identically zero since for (y,s) and (£, t) S <7((0, T)) we have Si(y, s;£, s) = 
S 2 (y,s;£„s) = 0. 

By (|3.5p we have 

(3.6) U(y,s;Z,r) = ^ < (A Ql - A Qa )r 2 (-, •; y, s), T{ (•, •; £, r) 

for all y,£ £ ft, < s < r < T. Denoting by Sl_D(t) := fi \ f e [0,T], 

we introduce a variation of the Hausdorff distance that, even though it is not a 
metric, we call it modified distance 

(3.7) d tl (t) = d tt (D 1 {t),D 2 (t)) 

= max < sup dist(ir, D 2 (t)), sup dist(x, Di(t)) > , 

I a;e9Di(t)naf2£,(t) xedD 2 (t)ndn D (t) j 



f € [0, T]. We point out here that trivially d ll {D 1 {t),D 2 {t)) < d n {Di{t),D 2 (t)). 
The following proposition shows the relation between the Hausdorff distance and 
dfj,, provided the priori assumptions (|2.7|l hold. We refer to [Al-DC, Proposition 
3.3] for the proof. 
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Proposition 3.2 Let D\{t) and D 2 (t) be two sets satisfying (|2.7p then for any 

te[o,T] 

(3.8) d n (dD 1 (t),dD 2 (t))<Cd ll (t) > 
where C depends on E and M only. 

We now give a proposition which connects the Hausdorff distance between the 
boundary of the inclusions and between the inclusions. 

Proposition 3.3 Let D\ and D 2 be two domains of class C 2 with constants 
E, po such that R ra \ Dj, j = 1, 2, is connected. There exists a positive constant 
C depending on E only such that 

(3.9) d H (D u D 2 )<Cd H (dD u dD 2 ). 
Proof. See Section IO 

Remark 3.4 By Propositions 13.21 and 13 .31 we have 

d n (Di(t),D 2 (t)) < Cd H (dD 1 (t),dD 2 (t)) < 

C'd^D^t)^^)) < C'd n (Di(t),D 2 (t)). 

Thus it turns out that the distances dn{Di{t), D 2 (t)), d-n(dDi(t),dD 2 (t)) and 
d^(Di(t), D 2 (t)) are equivalent. 

Proposition 3.5 Let J7i and VL 2 be bounded domains in K™ of class C 1 ' 13 with 
constants Rq, E and satisfying < MRq. There exist numbers d, ~p € (0, Rq] 
d ~p 

such that —— and depend on [3 and E only, such that if we have 
Ro Ro 

(3.10) d n {Th,Th) <d, 

then the following facts hold true 
i) Q,\ and Vl 2 are relative graphs and 

(3.11) 7(fii.^) <Cdn (Th,Th) , 
where C depends on [3 and E only, 

Hi) any connected component of £l\ C\Q. 2 has boundary of Lipschitz class with 
constants po, L, where po is as above and L > depends on E only. 



Proof. See [Vel Proposition 4.1.8]. □ 
A key ingredient for the proof of our stability theorem are fundamental 

solutions. We collect here some results we need. 

We shall denote by Tq(x — y,t — s) the standard fundamental solution of 

dt — A which is 

^-y,t-s) = ^ e-^, t>s. 
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We shall denote by T(x,t;y,s) the fundamental solution of the operator d t — 
div((l + (k — l)xg)Va;) (see |Arj). We recall that T satisfies the following prop- 
erties 

(3.12) T(x,t;y,s) = T*(y,s;x,t) V (x, t), (y, s) G Q, (x, t) ^ (y, s), 
where T* is the fundamental solution to — dt — div((l + (k — 1)xq)V x ), and 

(3.13) < r(x,t;y,a) < ^-^-^ e -^ X[j|+oo) (t), 

where C > 1 depends on k only. Furthermore we have also the following estimate 
for the gradient of T. 

Proposition 3.6 Let T(x,t;y, s) be the fundamental solution of the operator 
dt — div ((1 + (k — 1)xq)Vz). There exists C > 1, depending on k and E only 
such that 

C i — i 2 

(3-14) |V x r(a; 1 1; y, s)\ < _ e"&, 

(t-s) 2 

for almost every x, y £ M™ and t, s gi, t > s. 



Proof. See Section \£M □ 
In the sequel we need the fundamental solution of the operator £ + = dt — 
div((l + (k — 1)X+)V) where \+ = X{( a; ,t)6R"+ 1 -.x n >a}- We shall denote by 
r + such a fundamental solution. Also, we shall denote by r^j_ the fundamental 
solution of the adjoint operator of £+. Observe that T+(x,t;y, s) — T + (x 1 t — 
s] y, 0) and r^_(x, t; y, s) — T + (x, s — t; y, 0). Here and in the sequel, for a given 
function f(x , ,x n ), we shall denote by Tt> (f '(• , x n )) the Fourier transform of / 
with respect to the variable x' . Thus 

T c {f{-,x n ))= ( f(x',x n )e- ix '<'dx', 



for every (' G M n_1 . 

In [Is-Ki-Na] it has been proved some formulae for (T + (., x n , t; y)). The 
technique to prove such formulae is rather classical and lengthy. For this reason 
we display only the formulae that we need corresponding to the case in which 
x„ > 0, y n < 0. 

Case k > 1 . 

Denote by 



(3.15) E(C', x n , t; p) = cxp 



-t(k (k l)p)\C'\ 2 - \[^x n \C\J- p 



(3.16) F(C, y n ;p) = Im (A^p)^^^^) , 

where, for complex number z = a + ib, Im(z) denotes the imaginary part b of z, 
and 



(3.17) AM = ^ ' 



iy/k — It/1 — P + Vk^fp 
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Then 



(3.18) F c ,(T + (.,x n ,t;y,0)) = / \C\e- iy '<' E(C,x n , t; p)F((', y n ; p)dp, 



for every x n > 0, y n < 0. 
Case < k < 1 . 
Denote by 



G(C, y„, i; p) = exp [-t(l - (1 - %)|Cf + VT~fc2/„|C'IVP 
ff(C', x„; p) = Im (A 2 (p)e-""^^ c 'l 



where 



Then 



F c (r + (;x n ,t;y,0))= / \('\e-* y < G((' ,y n ,t; p)H(C,x n ; p)dp, 
Jo 

for every x„ > 0, y n < 0. 

Proposition 3.7 For every A S (0, 1] i/iere exist Ai, A2, A 3 S (0, Ao] such that 
for every h > i/ie following inequality holds true 



(3.19) 



\ dt \7 x T* + (x,t;-X 1 he n ,X 2 h 2 ) 
Jo Jr^ 

■ W x T (x, t; -X 3 he n , 0)dx 



> 



Ch nl 



where C, C > 1, depends on Ai, A2, A3 cmd /c only. 



Proof. See Section EH □ 
Through the paper we shall fix the value of Ai , A2 , A3 in such a way that (|3.19p 
is satisfied and we shall omit the dependence of various constants by Ai, A2, A3. 
In the following we shall often make use of this technical lemma whose proof 
can be found in [Rjl Lemma 3 pg. 15]. 



Lemma 3.8 Let a, f3 < ^ + 1 and a > 0. Then 
rt 



(i-T)- Q e~ 



d£dr 



where C depends on a, (3 and n only. 
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For t £ (0,T] fixed, we can assume, without loosing generality, that there 
exists O £ dDi(t) n 9fi£>(i) (for the sake of brevity we assume that O is the 
origin of W 1 ) such that 

(3.20) d^t) =dist(0,As(*))- 
Denote by 

p = min{d jLl (?),po}- 

Furthermore, denote by v(0,t) the exterior unit normal to dD\(t) in O pointing 
towards G(t). Now we introduce parameter S £ (0, 1] that we shall choose later 
on. We set 

(3.21) t 1 =t-\ 2 h 2 , y = \ 1 hv(0,t), yi = \ s hv{Q,t), 
where 

(3.22) < h < 5mm{p, Vt}. 

Notice that (I3.22| implies that t\ £ (0,7). By using (|2.7aj) it is simple to check 
that there exists C%, C\ > 1, depending on E only such that if 

(3.23) < 8 < ^ 
then, for every t £ [ti,t\, we have 

(3.24) dist(y, D x (*))>- min { \ x , A 2 , A 3 } h, 

(3.25) clist^i, Di(t)) > imin{Ai,A 2 ,A 3 }^. 



On the other side, using the inequality |Ve[ Proposition 4.1.6] 

(3.26) \dist(0, D 2 (t)) - dist( 0,D 2 (t))\ < —\t-t\, 

Po 

where Co depends on E and M only, for t £ [ti,t\ and by using the triangle 
inequality we have that there exists C 2 , C 2 > 1, depending on E and M only 
such that if 

(3.27) < 5 < ^ r 
then for t £ [ti,t] 

(3.28) dist(z, D 2 {t)) > ip, with z = y, y x . 



Proposition 3.9 Let {Di(t)} te R, {-D 2 (£)}teR oe two families of domains sat- 
isfying (|2.7|) and let Ai,A 2 ,A3 £ (0,1) be such that the inequality (I3.19|) is 
satisfied. Then there exist C , C > 1, and C , C > 1, C depending on k only and 
C depending on k, E, M, Ai, A 2 and A3 only such that 

(3-29) \u(y utl ;y,t)\>-^, 

for < h < i minjp, Vt}, where y±, ti,y,t, and p are defined in (|3.2ip . 
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Proof. See Section FO □ 



Theorem 3.10 (Two-spheres and one-cylinder inequality) Let A, A and 

R positive numbers with A £ (0, 1]. Let P be the parabolic operator 

(3.30) P = d t - di {j'dj) , 

where {a 10 (x, f)}™.- =1 is a symmetric nxn matrix. For £ € W 1 and (x , t) , (y , s) € 
R n+1 assume that 

n 

(3.31a) X\e< a * i ( a; .*)&&<A~ 1 |£| a 



</// 



d 



1/2 

(3.31b) ( ^ (o«(o: > t)-a«(v, a )) a ) < ^ (|x - y| 2 + \t - s\f 2 



-Let u be a function in H 2 ' 1 (Br x (0, -R 2 )) satisfying the inequality 

(3.32) |P«| < A + in ^ x (0 , 

Then there exist constants r\\ € (0, 1) and C G [1, +oo), depending on A, A and 
n only such that for every r%,r2i < r\ < r-i < rj\R we have 

(3.33) \\u(;R 2 )\\ LHBr2) < ^Hull^^o^jlluO.JZ 2 )!!^^), 



where 6\ = 



Proof. See Ne\ □ 



We can now start to prove our stability theorem. Before entering into details, 
we wish to warn the reader that, sometimes we use the previous auxiliary results 
(such as Lemma 13.81 or Proposition 13. 6p omitting some computations that are 
similar to the one contained in the proofs of Section 4. 

Proof of Theorem 12.71 We divide the proof of the theorem in two steps. In 
the first step we provide a rough stability estimate (see ()3.49j) below). In the 
second step we prove the logarithmic stability estimate (|2.10|) . 

Step 1. 

We shall denote by 

(3.34) E po ={x£ R" : p /2 < dist(af, Q) < p } 
and 

(3.35) fl P0 ={iel" : dist(x, ft) < p }. 
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Since U(y, s; £, r) is equal to when s > r (see Remark |3 . 1 [) . for (y, s) G S po x 
(0, T) we define 

(3.36) v(£, t) := U(y, s; £, r), (£, r) G 0((O, T)) := U t6(0) r)C(t) x {t}. 

Let /i be the one defined in (|3.22|> . with S < ~ and let Ai, A2, A3 G (0, 1) be such 
that inequality (|3.19[) is satisfied. Let x G 5 po be such that dist(x, M" \ S po ) > 
Po/8. Let us denote by 7 a simple connected arc in (fl w \ n.o(t)),Xi/i/2 = {x G 
f2 Po \ Q.o(i) : dist(x, 3(f2 Po \ 0£>(I))) > ^p}, connecting x to y, where y is 
defined in (l3~2T|l . By (|2~H]) we have 

Ce 

( 3 - 37 ) IK^)IU°°fs Po/2 (*)) < -jt, 



and by Lemma 13.81 and Proposition 13.61 we have 

C 

(3-38) ll w IU°=(e((o,T))) < — , 

Po 

where C depends on k only. It is easy to check that by (|2.6|) and (|2.7j) there 
exists C, C > C, depending on k, E, M such that for all x G (fi Po \ ^_d(<))ai/i/2 
and < h < ^ min{p, v^}, 



(•'-!• <7\,/, 1 : : B Xih (r) x (7-- J 



Cff((-oo,|. 



Since v solves the heat equation, we can apply Theorem 13. 101 along a chain of 
balls centered in points of 7. More precisely, let us define p = 77iAi/i/12, where 
771 G (0, 1) is defined in Theorem 13. 10[ and x%, i = 1, . . . , mh as follows: x\ = x, 
xi + i = j(ti), where U — max{i : |7(t) — Xi\ = 2p}, if |xj — y\ > 2p, otherwise 
let i — rrih and stop the process. We have m/, < CM , where C > 
is an absolute constant. By construction the balls Bp(xi) are pairwise disjoint 
and \xi-\-i — Xi\ — 2p for i — 1, . . . , m — 1 and \x m — x| < 2p. By an iterated 
application of the two-sphere and one-cylinder inequality (Theorem 13. 10p to v 
with R = X\h/A, T2 — 3p, r\ — p over the chain of balls £> p (xi), since we have 
B ri {xi+\) C B r2 (xi), i = 1, . .., m/j, by (|3.37p and (|3 . 38[1 we have 

(3-40) (±- [ v 2 (U)dt 




\'2 JB r2 (y) 

where S2 G (0, 1) is an absolute constant and C depends on a priori data only. 
From now on, in order to simplify the writing and since the case e > 1 is trivial, 
we shall assume that e G (0, 1). By standard regularity estimates [Li] and taking 
into account (I3.37|) and ()3.38|) we have 

(3-41) l|V«(-,t)|U-( flra/a Cy)) < 

where C depends on k only. Recalling now the interpolation inequality (see 
|Al-Be-Ro-Vel (5.29)]) 



(3.42) \\g\\ L oo (Br) < C [\\g\\ L oo (Br) + r||Vfl|U- (flp) ] "+ 2 (r~" g 2 
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for every r > 0, where C is an absolute constant, by (|3.41[) and (|3.40p we have 

(3-43) \\v(;t)\\L^Br 2m) <^ h :=^\e), 

where C depends on the a priori data only. Now defining w(y, s) = U(y, s;y,t) 
and taking into account (|3.43p we have 

{ d s w(y,s)+A v w(y,s) =0, inS((0,T)), 

{ \w(y,t)\<wg\e), (2/,s)€S P0 x (0,T). 

Now we want to estimate from above In order to obtain such an 

estimate we argue as before, but here, instead of ()3.38p . we use the inequality 

C 

( 3 - 44 ) \\ w \\L~(gW([ tl ,t\)) < ^ 

where G W ([h,t\) = {(x,t) E 0((O,T)) : h<t<t, dist(x, fl D (t)) > S 2 h}, S 2 = 
g min{Ai, A2, A3} and C depends on k only. Inequality (I3.44| is a consequence 
of Proposition 13.61 and Lemma 13.81 Notice that by virtue of (13.241) and (|3.25|) 
we have (yijii), (y,t) € Q^ h \[t\,t\). Therefore we obtain 

(3.45) \U( yi ,h;y,t)\ = \w( yi ,h)\ < —e^ h , 

where s 2 € (0, 1) is an absolute constant and C depends on the a priori data 

only. Now we introduce some notation. We set 5(t) = min l| and h/po = q- 

By Proposition 13. 91 and (|3.45j) we have that there exists C5, C5 > 1, depending 
on k, E and M only such that 

1 _:_ K(*) 



(3.46) 1 < C 5 £ S3 , for every q < — min < - J ^- L ,5(t) } , 

C5 I Po J 

where S3, S3 G (0, 1), depends on M only. We distinguish two cases 

i) (?) > min{ Vt, p } ii) (?) < min{ a/?, Po } ■ 

If case i) occurs we choose q such that s| = | loge| -1 / 2 , that is 

1 

I log s 3 1 



q = Qe ■ = 



log|loge|-V2| 



Denote by £*(?) the least upper bound of the set {e € (0, 1) : q e < S(t)}. By 
(|3.46p we have 

l<C 5 exp{-|loge| 1/2 }, 

which, for < e < e**(i) := min |e* (?) ,e^' logC5 ' |, yields to a contradiction. 
Thus, if < e < £**(?), case i) cannot occur. 



15 



Let us consider now case ii), that is d^{t) < min |vl, po|. By (3.46) we have 

l<C" 5 cxp|-s| |loge||, 



GZHn (t) 

for every q < — . Now, if 



Pa 



(3.47) 



d M (*) < 2C 5 p 



log |loge| 



then we are done. On the other side, if 

(3.48) (i) > 2C 5 p a 

let us denote by 



log | log e | 



log I | logs 

and by e (*) the least upper bound of the set {e € (0, (I)) : g e < £ m}. Now, 
for < e < £q (i) , we choose q = q e and by (|3.46p we have 



l<C 5 exp{-|log £ | 1/2 } 



Since the last inequality yields to a contradiction whenever < e < £o (t), we 
have that if < e < £o (t) then (|3.48p cannot occur, so inequality (|3.47|) holds 
true . Finally, by using Proposition 13.21 and 13.31 and Remark 13.41 we have 



(3.49) d n (D 1 (t),D 2 (t)) <2C 5Po 
for < e < e (t). 



log |loge| 



1°8 S3 



-Xj n 



Step 2. 

In order to prove estimate (|2.9p we apply Proposition ^. 51 to := 0\ Di(t), 
i = 1,2, and Rq = po- Indeed by ()3.49p we have that, for e small enough, fli 
and f^2 are relative graphs. More precisely if < e < mm{eo(t), d} (where d is 
defined in Proposition 13. 5p then there exists r$ > (ro/po < 1 depending on E 
only) such that 

(3.50) B ro (0) n Aft) = {x e Bro(O) : x n > Vl {x')}, i = 1,2, 

and ||(/3i — c^2 1| (s; ) < Cpoc(£); where C depends on only. By (|2.7ap and an 
interpolation inequality |Al-Be-Ro-Ve|. (5.30)] we have that \\fi — f2\\c 1 (B' rg ) < 

Cp {a{s))P, ]3 e (0, 1). Thus, with eventually a rigid transform of coordinates, 
provided we pick a smaller tq, there exists £o > 0, depending on E only, such 
that for e < So we can assume that |Vy>i(0)| = 0. In the sequel we continue to 
denote by £o(i) the number min{eo(^), £o}- Let us define, for a unit vector £ and 
< a < tt/2 



(x — z) • C 

x G B r Jz) : — : : — > cos a 

\x-z\ 
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By (|3.50p we have that C(0, v, a, ro) C G(t), where a, a G (0, 7r/2), depends 
on E only. Let us denote ~p — 1+ g° n - , <5* = 1 an( i Po = 7^ cos § • We have 

that S((0,t),is, j,5* ,~p) C <?((— oo,t]), where we set v — v(0,T), for the sake of 
brevity, and 



= G M n+1 : z G C(0,v,^,p), t~S*(x ■ vf < t < <| 

and ^ = u(0,t). We want to estimate v(y,t) = U{y, s;y,t) when (y, s) G S po x 
(0, T), s < i, where v solves 

' d T v-Asv = Q, for (£,r) GS((0,T)), 

(3.51) J «(e,r)| T < a = 0, 

k IMU°°(e po x(o,t)) < e, 
where £ = Ce/pQ, where C depends on k only. Let us fix (y, s) such that 

(3.52) (y, s) G S po x (0, T) = { l£ K" : ~p < d(z, n) < 7 -p } x (0, T). 

By Proposition (|3.6p and Lemma (|3.8p . taking into account the last relation of 
(|33Tj) and by (j3"32"j) we have 

C 

(3.53) IMIl»(s(o,t)) - 7^ + e:=H. 

In the sequel we continue to denote by v the trivial extension of v. Notice that, 
by (|3.53|) . we have 

\\ V \\L^(S({0,t),v,a/2,5' ,p) ^ H. 

Denote by 

ot\ — arcsin ^min |sin — , S* ^1 — sin — ^ |j , 

Ml = ^— — : ,wi = z + HiC,Pi = -tViVi Sln «i. 

1 + sinai 4 

where rji G (0, 1) is defined in Theorem 13. 101 We have 

(3.54) dist (wi, dQ (I)) > min {po — |it?i — z\ , |iui — z| sin a} = porj, 
where 

f cosf 1 

?7 = min < 1 — 



1 + sin a 1 + sin ai 1 + sin a\ 1 + sin a 

Now (I)) ppj? is connected and, by (|3.54p . w\ G (i))pa5- Therefore by an 
iterated application of the two-sphere and one-cylinder inequality (see also [Vel 
Proposition 4.1.1]) we get 

(3.55) ( pr" / v 2 (U)dti) KCe^H 1 -*, 
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where S4, S4 G (0, 1), is an absolute constant and C depends on the a priori data 
only. Denote 

Hk = a k ~ 1 n\, w k = fj, k p p k = a k ^pi 
dk = fJ-k- Pk = a fe_1 /il(l - sinai), 

where a = \ }'" smai _ p or every k > 1, the following inclusions hold true 

1+^771 smai J — ' 

(3.56) B Pk+1 (wk+i) C B 3pk (w k ) C B 4r -i pk (w k ) C C(0, i/,oti, r ) 



(3.57) B 4 „-.(%) x (i-(4^V) 2 ,*] C5((0j),i/,a,<j*,r ). 



Let us consider h defined in (|3.22|) . We further assume X\h £ (0, d±]. Let k be 
the smallest positive integer such that dk < \\h. We have 

(3.58 n : <k-l< : hi. 

|loga| |loga| 



Denote ^ ^ 

\ J B Pj { Wj ) ) 

By the Theorem EUDl (EHJ), (l3~5?l) and since 



j = 1, ...,fe. 



V J B3 Pj ( Wj ) ) 

we obtain 

(3.59) a^KdH^-^af*, j = l,...,k-l, 
where 0* = Co i og (4 /m) ■ By iterating (|3.59|) we get 

(3.60) a l<C^H 2 
By Lemma T3. 81 and Proposition 13. 61 we have 

(3-61) IM|L~(Q Pr ) < 

where = B^m^) x (f — p|,t]. By standard regularity estimate and (|3.6ip 
we get 

( 3 - 62 ) llv«(- s t)|| £ » (B (wr)) < JL. 

k PkPo 

Now by using interpolation inequality (|3.42p . (I3.55[) (|3.60[) and (|3.62[) . taking 
into account that y € B p _/ 2 (. w k)^ we have 

(3-63) \ v{ y- t) l<c(^Y 2 H ((± 



Pk 



H 
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where 9* € (0, 1) and s 5 e (0, 1) (absolute constant) and C depends on the a 
priori data only. Now evaluating k in terms of h and recalling that w(y, s) — 
U(y, s;yj) = v(y,t) we have, for every (y, s) € S po x (0,T) 

I log e, | 

(3.64) Hl£- ( i W)XC0 , T) j<^(^) ' 

where C depends on the a priori data only. Arguing as above to estimate 
\w{y\,ti)\ and recalling that w(yi,t\) =U(yi,ti;y,t) we have 

(3.65) \U( yi ,ti;y,t)\ < ^e^) , 

where C and B depend on the a priori data only. Finally, using Proposition l3.9l 
and proceeding as in Step 1 we obtain (|2.10jl . □ 



4 Proof of the Auxiliary Results 
4.1 Proof of Proposition EOI 

Proof of Proposition 13.31 We recall that, for a given a subset A of R n , we 
denote by [A] e = {i£l" : dist(x, A) < e} , (A) e = {x e A : dist(x, dA) > e) 
and [dA] e MP : dist(x,&4) < e}. We remark that [A] e \ (A) e = [9A] S 

and (A) E C Ac [A] E . 

Let d := d H (D 1 ,D 2 ) and r := d H (dD 1 ,dD 2 ). If d = then fTS]) holds trivially. 
Assume <i > 0. Without loss of generality, we can assume that there exists 
x S D\ such that d = dist(x, D 2 ). Since e? > we have that x £ D 2 and 
therefore d = dist(x, dD 2 ). If x € dD\, then (|3.9p is trivially true. Assume 
x € -Di \ D%- We have for every x € D\\D 2 

dist(x,<9D 2 ) = dist(x,Z> 2 ) < dist(x,Z> 2 ) = dist(x, dD 2 ). 

Thus for every x G Z?i \ £) 2 we have 

dist(x,&D 2 ) < dist(x,&D 2 ), 

that is x is a maximum point in the set D\ \ D 2 for the function dist(-, dD 2 ). 
In the set A = Int ([dD 2 ] PB / E ^J \ dD 2 , the function dist(-, dD 2 ) is C 2 and 

(4.1) |V K dist(x,aL» 2 )| > VxeA. 
Since x is a maximum point and x ^ 9£> 2 , by (14. ip we have 

(4.2) dist(x,aD 2 ) > ^. 

Otherwise, recalling that x is a maximum point of dist(x, dD 2 ) interior to D\ \ 
D 2 , if dist(x, dD 2 ) < po/E wc should have V K dist(a;, dD 2 ) = contradicting 
(|4.ip . First let us assume r be such that 

(4.3) r<mm |_ _j. 
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We can write R™ = (D 2 ) r U [dD 2 } r U (R™ \ [D 2 ] r ), By (g^) and flOJ) we have 
x ^ [9Da]r- Since (Z?2)r C -D2 and x ^ D2 we have that x G R n \[D 2 ]r- Recalling 
that r < po/2 and R™ \ D2 is connected, we have that R" \ [D 2 ] r is connected. 
Thus there exists a continuous path 

7 : [0, 1) R™ 

such that 

(4.4a) 7 ([0, 1)) c R" \ [D 2 ] r , 



(4.4b) 7(0) = x lim 7 (t) = 00. 

t->i~ 

Since d H (dD 1 ,dD 2 ) = r and dL>i C [<9L> 2 ]r C [D 2 ] r , by (H^aj) we have 

7([0,l))nflDi = 

which is a contradiction since 3: G Z?i and Z?i is bounded. Thus we cannot 
connect x and 00 with a path that does not intersect dD\. Hence x G dD\. 
Thus if du(dDi,dD2) < 5poi with <5 = min{-g,i}, (|4.3[) is satisfied and we 
have 

d H (p 1 ,D 2 )=^st(x,dD i ) < r. 
On the other side, if d-n(dDi,dD 2 ) > 5po we have trivially 

(4.5) d n (D, ,D 2 )< 2diam(tt) < 2dl ^ m (") ^ {dDl _ dDi) 

and the proposition is proven. □ 

4.2 Proof of Proposition 13.61 and Asymptotic Estimates 
for the Fundamental Solution 

We shall make use of the following regularity theorem, whose proof can be found 
in [La Ri-Urj . [La-So-Url Ch. Ill, Sec. 13]. 

Theorem 4.1 Let X, M and r be positive numbers with X £ (0,1]. Let u G 
H x ^{B' r x (— r, r) x (— r 2 ,r 2 )) be solution to 

(4.6) div (A(x, t)V x u) + b(x, t) ■ V x u ~ d t u = 0, 

where A(x, t) and b(x, t) are respectively a symmetric nx n matrix and a vector 
valued function satisfying the following conditions 

(4.7a) A|£| 2 < A(x,t)£ ■ £ < A _1 |^| 2 , 

for all (x,t) G B' r x (— r,r) x (-r 2 ,r 2 ) and for all £ G R", 

n n 

(4.7b) r ^ ||9jA||i«,( B ^ x (_ r) o) x (_ r 2 )r 2)) + r ^ ||9iA||ioo( fl / x (o,r)x(-r 2 ,r 2 )) 
i=l t=l 

r2 ||^t^lU°=(S;x(-r,r)x(-r 2 ,r 2 )) < -W, 
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(4.7c) r||6|| L oo (B^x(-r,r)x(-r 2 ,r 2 )) < M. 

Then there exist positive constants (3 € (0, 1) and C such that for every p < § 
and all (x,t) e B' r _ 2p x (— (r — 2p), (r — 2p)) x (— r 2 + 4p 2 ,r 2 ) </ie following 
inequality holds 

(4.8) pllV^ullioo^^xC-p+^p+z^xC-p^+t,*)) 

+p' 3+1 [V a; lt] / 3 ; ( B / ) ( :I ./)x(-p+2:„,p+a;„)x(-p 2 + tt))n(S^x(-r,0)x(-r 2 ,r 2 )) 
+P /3+1 [Va; M ]/3;(B^(a;')x(-p+2:„,p+a;„)x(-p 2 + tt))n(S;x(0,r)x(-r 2 ,r 2 )) 

< J / u 2 ^,T)d^dr\ . 

Here (3 depends on n only and C depends on X, M and n only. 



Before proving Proposition ^. 61 we give the following estimate which is needed 
in the proof. We recall that Q p (xo,to) = B r (xo) x (to — p 2 ,to)- 

Proposition 4.2 There exist constant C > 1 and < Si < 1 depending on k 
and n only such that the following inequality holds. 

(4.9) / \T(x,t;Z,T)\ 2 dxdt < C P " - e~W^, 

jQ P (x ,t ) (*0-T)« 

w/iere p = ft [|x - £| 2 + <o - r] 1/2 . 

Proof. From the inequality ()3.13|) we have 
(4.10) 



/ 

Jo 



,T{x,t;i,T)\ 2 dxdt<C [ 1 e ^I^T X[t +oo) dx dt, 

(x ,to) JQp(xoJo) \ l T ) 



where C± depends on k and n only. In what follows we denote by / the integral 
at the right-hand side of (|4.10p . We distinguish two cases 

i) t - p 2 < t < t , 

ii) r < t - p 2 . 

Let us consider case i). It is easy to see that there exists an absolute constant 
C > 1 such that 

(4.11) C _1 p < \x - C| < Cp VieB,,(io). 
By (|4~TT|) we have 

(4.12) I < c nP n f ° S - n e~^ds, 

Jo 

where c n is an absolute constant depending on n only and C2 depends on k and 
n only. Now if < to — r < ;^tj, being s — > s ™e C2S an increasing function in 
(0,^), by 612) we get 

(4.13) I< ^ e-^T. 

(t - r) n 1 
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Otherwise, if -^r < t — r < p 2 then since 



max 

(0,+oo) ' 



and now t — r is of the same order of p 2 we have 



p- n I < d 



l e ^ds < 



C 



(t - t)«-i 



e C 2 (t — r) _ 



By the last inequality and (|4. 1 3(1 we get the Proposition in case i). 
Let us consider now case ii). It is easy to see that 

(4.14) 6p 2 < \x-£\ 2 +t-T < 60p 2 , 

for every (x,t) € Qp(xo,to). Moreover, denoting 

i_ 

J e 

ip — 



e c i('- t ) 

M„ = max < — r— : (x,t) £ Q p {x ,t ) 



(t- T y 

and taking into account (|4.14p we get 
(4.15) 



where C depends on n only. Now, since r < to — p 2 we have 

\xo-t\ 2 



(4.16) 



i - t 



< 4. 



Therefore by (|4.15[) and (|4.16p we get the Proposition in case ii) as well. □ 



Proof of Proposition 13.61 Let U be a solution of the equation CU = 0, 
where C = d t — div(l + (k — 1)xqV). We recall the following regularity estimate 
(see |La-So-Urp 



(4.17) 



c 

n + 4 
f 2 



\ 1/2 

_ U 2 (x,t)dxdt ) 

Q2r(^,t) / 



where Q+(x,t) = Q r (x,t) (~1 Q (we recall Q = f(K)) and Q r (x,t) — Q r (x,t) \ 
Q + (x,t). Applying (|4.17|) to the function P(-, -;£,,t) we get 



(4.18) ||vr( v ;£,r)IU (Q ± (xoito)) <4L 

P 2 



-I 1/2 



\T(x,t;i,T)\ 2 dxdt 



1 



9 = 4 [l^o-€| 2 +io-^] 



1/2 



where 
(4.19) 

Applying Proposition 14. 21 to the right hand side of (|4.18|) we have 



|Vr( v ;£,T) 



\L°°(Qf{x ,t a )) 



< 



c 



(to - r) 



n-l 



e c{t -T) 



1/2 
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Since 



1 1 

< 



16p 2 ~ t -T 



we obtain filty . □ 

In order to state the next theorem we introduce some notations. Let ip : 
B' x (— PqiPo) — * ^ suc h that it is difFerentiable with respect to t and Xi, 
i = 1, . . . ,n — 1, it is twice difFerentiable with respect to Xi, i = 1, . . . , n — 1, 
and df(p is difFerentiable with respect to Xi, i = 1, . . . , n — 1. We assume that 



(4.20) 
and 



^(0,0) = |V a ^(0,0)| = 



(4.21) p Q \\D x ,(p\\ L ao( B > x (_ /5 2 )p 2)) +Po||9t¥>||ioo( B / X (-p 2 ,p 2 )) 



We shall denote by 

Qp, P0 ={xe B po x (-po»Po) : x « > 

and by Ta+ (x, i; y, s) the fundamental solution of the operator dt — div((l + 
(k — l)xn+ V), that is 

dtT Q+, P0 Mil/,*) 

- div ((1 + {k - 1)Xq+ pQ ) vr Q+ P0 ( x : *5 s i) = ~ S ( X ~ V, 1 - s )-, 

where (y,s) £ R" +1 . 

Theorem 4.3 (Asymptotic Estimate) Let ip and Tq+ (x,t;y,s) as above. 
Then there exists a constant C > 1 depending on n and E only such that 



(4.22) T d+ (x,t;y,0)-T + (x,t;y,0) 



y\ 2 +tY' 2 e-^cf 



Po 



t n/2 



(4.23) V i r + pQ (a;, i; y, 0) - V x r+(a;, i; y, 0) 



Po 



t n/2 ' 



where [3 is the one defined in Theorem \4 ■ 1\ depending on n only, for all 
(x, t) G Q Vi m n | (a:, t) e R" +1 : t > 0, x n > ^(l^'l 2 + *) 



and y = y n e„, y n € (-p o /C,0). 
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Remark 4.4 Theorem l4.3l provides an asymptotic estimate for the fundamental 
solution T(x, t; y, s) when [x, t) and (y, s) stay on opposite sides of the interface 
(given by the graphic x n — ip(x',t)). Our crucial requirement is that (y,s) 
approaches the interface in a nontangential way. 

Proof, of Theorem Let 9 be a C°° function on K such that < d < 1, 
9{s) = 0, for every s G E\ (-2,2), 9(s) = 1 for every s G (-1, 1) and |0'(s)| < 2 
for every sel. 

We define new variables by (£,t) = ^(x,t), where ^(Xjt) = (&(x,t),t) and 

r r = x', 

( T = t, 

where n = p min{i, 3^}. 

Sometimes, for the sake of brevity for a fixed t G (— PgiPo) we denote by 
<£>(*)(•) the map $(-,£) and by the graph of tp(-,t). It is not difficult to 
check that and $ have the same regularity properties of ip and they are 
diffeomorphisms (that preserve orientation) of and M. n respectively. We 

denote by $~ 1 (-,f) the inverse of $ < - t - ) (-). The following properties hold: 

(4.24a) $W n (B^ x (-n,n))) = {a: G X (-n.n) : x„ = 0}, 

(4.24b) *(x,<) = 

V(x,i) G M" +1 \ ((^ ri x (-2ri,2n)) x (-2r?,2r?)) , 

(4.24c) C-^asi-xal < |$ (t) (^i) " * (t) (^)| <C|a?i-x a |, Vxi,aj a eM n , 

(4.24d) |$ (t) (ir) - x| < — Id 2 , Viel", 
Po 

(4.24e) |£> x $W(x)-/| < — \x\, VieR", 

Po 

where C, C > 1, depends on only, / denotes the identity matrix and D x ^^ 
is the jacobian matrix with respect to variable x. For y n G (— 9-,0) and a G 
(— r^,r^) we denote y = y^en and 77 = $(j/,<t). Furthermore we shall use 
the following notation f (£, r; 77, cr) = r(* _1 (£, r); cr)), and 7(£,r) = 

det J(£,t), where J(£,r) = (£> X $)(*~ 1 (C, r)). We have that f(£,r; 77, er) is a 
solution to 

(4.25) div r)V 4 f ) + C(£ r)V s f - e r f = - 7 fa, <t)i$(£ - ry, r - a), 

where B(£,r) = (1 + (k - l) X +)S(?,r), B(£,r) = (J(£, r)) (J(£, r))* and 

^ I r) = J^r)^^-fg^V a (C,r). 

Since we want to study the asymptotic behaviour of r(£,r;?7,0), we shall 
denote f(£,r;r7,0) byf^r;^). 
By (gHU), we have that 

(4.26) 5(0,0)=/ and \\B\\ Laa(ax( p tT)) + Po[£]i,nx(o,r) < C 
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where C depends on E only. Denote by 

(4.27) R{i, r; V ) = f (£, r; ry) - 7(7?, 0)r+fo r; r?), 

where r + (^,T;?7) = r + (£, t; 77, 0) is the fundamental solution to the operator 
div ((1 + (k - l)x + )V 5 ) - d T . We have 

div ((1 + (k- l)x + )V e i?) - d T R = r; ??), 

where 

F&tjtj) = -C(£, r)V € f (£, r; tj) 

+div ((1 + (fc - 1) X +)(J - r))V f f (£, r; 7?)) , 

notice that F(£, t; 77) = if r < 0, 

(4.28) #(£,7-577) = 0, for r < 0. 
Therefore, lArl 



R(^t; V )= / Ffcaj^r+for; if r > 0. 



We have 

(4.29) 

where 

(4.30a) 
and 

(4.30b) J 2 = 



Ji = 



JB 2r 



\R{^r-vi)\<J 1 + J 2 



C(C, s) V c r(C, s; 77)r+ (C, r; C, 



JB 2r 



(l + (k-l) X +) (I-B(C,s)) 



V c f(C,s;77) • V c r + (£,T;C,a)dC<fa 



By Proposition [32] and Lemma T3. 81 we have 



(4.31) 



Ji < 



-e c i T . 



Pot" 



where C, C%, C > 1, C\ > 1, depend on only. By Proposition 13.61 and (|4.26|) 
we get 



(4.32) 
with 

(4.33) 



Ji < </2,l + ^2,2, 



err n ic-^i 2 . , 2+1 ie-ci 2 

J 2 ,i = — / / s"e c 13 (r-s) ~e c i^-^d(ds, 
Po Jo JM n 
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and 
(4.34) 



Cf T f a±i \<-v\ 2 , N n+i ie-ci 2 

•ha = —\ / |C|s ~e c i s (r-s) ="e c i<-—°> d(ds, 

PO Jo JR n 



where C,C\, C > 1, Ci > 1, depend on E only. 
By Lemma 13.81 we obtain 



(4.35) 



Ja,i < C- 



_2_1_ i 

T 2 "t" 2 || - t; I - 



-e c ^ 



Let us consider now Ji,i- Performing a change of variables we get 



Ji 



2,2 



C_ (VCix e~^T 
Po V 2 



/2 



< 



PO T 



n/2 



^( T (l-A)A) 1 /^ + A (e-7 ? )+^ 



m 2 + T) 1/2 + \a-v\ 



dzdX 



where C depends on _E only. Now, denoting by C2 = max s 1 / 2 e 2C "i s , we 

se(o,+oo) 

have 



e c i- 
r™/ 2 



1 /:c-r/ 2 ^ 



C* 2 



Thus 



(4.36) 



/7 p 2C'it 

^ 2 <-^r7^(ICI 2 



PO T 



l/2 



Now since rj = e n r] n , rj n < and£„ > 0, w e have \£ -r]\ 2 = |£| 2 - 2?7„£„ + |?7| 2 > 
|£| 2 . Such an inequality and ([428| . (|43T|) . ([435]) . (|436l> give 



e c 



(4.37) TOr;»7)|<£x^ xp ,. +00) T „ /2 

for every £ £ and r G (0, 4r 2 ), where C, C > 1, depends on _E only. Let <5i 
be the constant defined in Proposition (Si £ (0, 1)) and, for fixed £ <G -B^/ 8 > 
??n € (-ri/8,0), r] = e n r) n , f G (0, (ri/8) 2 ) denote by 



h= j[\I-v\ 2 +r} 1/2 - 



We have 



div (%, r)V € f ) + C(t r)V ? f - d T f 



0, 



in B^ /2 (£ ) x (£ n - /i/2,£„ + /i/2) x (t- (h/2) 2 ,r], where 5 and C are defined 
above. Therefore by Theorem 14. f I and Proposition 14. 21 we get 



(4.38) 



CI 
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where where Q = B' h u{£) x + h/4). Since a similar inequality holds true 

for V ? r + (-,T;?7), by (|Q7|) we obtain 

(4.39) W,?;^^^-^. 

In (|4.38p and (|4.39p . C, C > 1, depends on _E only. Now we recall the following 
interpolation inequality 

(4.40) ||V/|U~ (Q) < C (||/||£? (Q) |V/|^ + i||/|U~ (Q )) • 
Since (|4.37jl easily yields 



(4-41) \\R(;r;r,)\\ L ^ Q) < —t-Z-h, 

PO T n/Z 

where C, C > 1, depends on E only, we obtain by (|4.39[) and (|4.40p 



C ( h 



-1+7 



(4.42) |V^(^r;7?)| < - - -^/a". 

for every J G S+ /gJ rj = e n rj n , rj n G (-n/8,0), r G (0, (ri/8) 2 ], where C, C > 1, 
depends on E only. 

Let us go back to the original coordinates (x,t). First of all let us estimate 
the function g defined by 

(4.43) g(x, t; y) := (x), i; §<°) (y)) = (a?),*; e„y„). 

To carry out the estimates, up to the end of the proof, we always consider x 
and y n such that x G Bs Po (Spoe n ), y n G (— (JpcbO), where J, <5 G (0,1), may 
change from line to line, but it shall depend on E only. Notice that for every 
x G Bs Po (5poe n ) we have > 0. Also notice that 

(4.44) \x\ < \x - e n y n \, x G B$ po (5p e n ), y n G (Sp ,0). 
By such an inequality and (14.24dl) we have 

(4.45) \^*Hx) - x\ < —\x - eny n \ a 

Po 

for x G Bs Po (Spoe n ), y n G (— <5po,0) where C depends on E only. By (I4.45[) . 
(|4.44|) and the triangle inequality we have 

(4.46) C~ x \x - e n y n \ < |$ (t) (x) - e n y n \ < C\x - e n y n \, 

for x G Bs Po (5poe n ), y n G (— £po,0) where C, C > 1, depends on E only. By 
(IPil . (f4~4% . lf£23e)l . (14^31) . (14441) . d446)l we obtain 

I a - e„j/„| 2 + tl 17 ~ 



e ct 

(4.47) \g(x,t,y)\<C 



Pi 
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and 
(4.48) 



\V x g(x,t;y)\ < 



C e-^^ 1 n\ x - en y n \2 + t] i/2 

Po t^{\x-e n y n \2 + t}^\ po 



1+0 



for every x £ Bs Po (8poe n ), y n G (—Spo, 0) where C, C > 1, depends on i? only. 
Recalling the definition of g we have that 

(4.49) 7 ($ (0) (y)) (rfo t; y, o) - r+(s, t ; y , o)) = 

9(x,t;y)- (l-7(* (0) (2/))) r + (a;,t;y,0) 

- 7 (V 0) Q/)) (r+Car.tjv.Oj-r+^^a:),*;!/^)) . 

Now for x G Bs Po {Sp e n ), y n £ (Sp o ,0) we have 

(4.50) \y\ < \x - e n y n \, 

so such an inequality, (|4.26|) . (|4.24cp and (|3.13[) give 

(4.51) |(l- 7 (^(y)))r + (x ) *;y,0)|<C^- &nVn? " ' 



p8 



t™/ 2 



for every x £ Bs Po (SpQe n ), y n £ (—6po,0) where C, C > 1, depends on only. 
In order to estimate from above the third term at the right hand side of (|4.49|) 
we use the mean value theorem. By such a theorem, (|4.26l) and Proposition 
we get 

(4.52) | 7 ($ (0) (2/)) (r + (*W(x),i i y 1 0)) -T + (x,t;y,0) 

l~-<=„y„| 2 



n + l > 
t 2 



where x — x + A($(*)(x) — x) for a suitable A £ (0, 1) and C, C > 1, depends 
on .E only. Now, by triangle inequality, (|4.24d[) . (|4.44|) we have 

(4.53) \x - e n y n \ > \x - e n y n \ - \x ~ x\ 

>\x- e n y n \ (l - — \x\ ) > \\x - e n y n \, 
V Po J 2 

for every x £ Bs po (5pae n ), 6 small enough and depending on E only. By in- 
equality jUMJ, (|g^E5jl . (fOT]) . (j43Tj) and ((43211 we obtain 

(4.54) |r+(x,i; 2 /,0)-r + (<i>( t )(x),i; 2 /,0)| 

N-e„t/„| 2 +t)\ 1/2 



< C 



e 



for every x G B Spo {5p Q e n ), y n £ (Sp o ,0), t £ (0, (ri/8) 2 ), where C, C > 1, 5, 
< 5 < 1, depend on E only. 
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We finally esti mate \ \7 x T(x, t; y,0) - \7 x T + (x, t; y, 0)| for x £ Bs Po {5p e n ), 
y n € (-5/5o,0). By g^SJ), (|OU|) , gUJ), and Proposition [32] we have 



(4.55) | V.rfo «; y, 0) - V x T + (x, t; y, 0)| 

Ce"^ /[|z_ en y n |2 + t] i/2 



< 



Po i" /2 V Po 



+ 1 



C 



.(x,t;y,0)-V x {T + (^(x),t;y,0) 



for every x € B$ Po (5p e n ), y n € (-<fy> ,0), i € (0, (ri/8) 2 ), where C, C > 1, 
depends on £7 only. Let us consider now the last term at the right hand side of 
(@~55|). We have 

(4.56) | V x T+(x, t; y, 0) - V x (l+($« (x), i; y, 0)) | < H^x, t; y) + H 2 (x, t; y), 
where 

H x {x, t- y) = C\V x T + (x, t; y, 0)| |7 - A^^x)] 

and 

H 2 (x, t; y) = C\D x ^\x)\ \ (V X T+) ($(*>(*), t; y, 0) - V x T + (x, t; y, 0)|, 
where C depends on E only. By Proposition ^. 6[ (14.24c[) and (|4.44p we have 

_ \ x-enVn I 2 _ I ^ — ^nVn I 2 

e e 2ct 

(4.57) H 1 (as,t,»)<C|a;-e n y n | ^ < C- 



t z±i t n,2 > 

for every x € B Spo (Sp e n ), y n E (-6p ,0), t e (0, (ri/8) 2 ), where C, C , depend 
on E only. To estimate from above the function H 2 we apply Theorem 14. II Let 
Si be defined as in Proposition ^. 21 and let us denote 



p=^[\x-e n y n \" + t}^. 

By (|4.24dp and (|4.44p we have that there exists S € (0, 1), depending on E only 
such that 

(4.58) |$^(x)-x| <^p, 

for every x G Bs Po (5poe n ). Now Theorem 14. II and Proposition 14. 21 yield 

(4.59) [V a r + (.,t; 1 ,,0)] / ,. B / (a / )x(a! o taB+/ ,) <Cp-( 2 +«- ^ , 

t 2 

where x„ = max{0, x° — p} and C depends on E only, x € Bs PQ (5poe n ), y n £ 
(-5p ,0), t e (0, (5po) 2 )- By (14.24d|) . (j4~59l) and (1435]) we obtain 

^ _ \ x ~ e nVn I 2 

(4.60) ^2(x,t;iy)<-|x|V( 2+ «- e 



Po t~ 
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for every x £ B Spo (5p e n ), y n £ {-5p 0) t £ (0, (<5p ) 2 ), where C, C > 1, 
depends on £ only. Finally (T04")) and (|^THU|) yield 



Po W 2 V Po 
The last inequality, (037]), (j4~55]) . fl4~56"l) give 

\V x T(x,t;y,0)- V x T+{x,t;y,0)\ 



< 



Po t n / 2 \ po 



for every x £ Bs Po (5poe n ), y n £ (—5po,0), t £ (0, (8po) 2 ), where C depends on 
E only. □ 

4.3 Proof of Proposition 13.91 

Proof of Proposition [3771 First of all, let us observe that 

(4.61) r (£, r; -A 3 e„, 0) = h n T {h^ h 2 r; -X 3 he n ,0) 
and 

(4.62) r;(£, r; -A ie „, A 2 ) = h n T* + (h£, h 2 r; -Aifte n , A,/! 2 ). 

Indeed (|4.6ip is a trivial consequence of the definition of To Concerning (I4.62[) . 
it can be proved as follows. Denote by y^ = —Xihe n , = A 2 /i 2 we have 

(4.63) / (r* + d t ip+(l + (k-l) X +(x))V x T* + -V x <p)dxdt = <p(y( h \sW) ) 



for every ip £ C£°(R n+1 ). In T* + denotes the function T*(x, t; yW, sW). 

Now for an arbitrary function ^ € C£°(R n+1 ) put p(x, t) := if) (f , in |(C53) . 
In the obtained integral we perform the change of variables x — h^, t = h 2 r. 
Thus, taking into account that x+{hQ — X+(0> h > 0, and denoting by Wh{£, t) 
the right-hand side of (|4.62|) we have 



{Whdri, + (1 + (k - l)x+(0)V£W fc • V £ V) ^dr = ^(-A ie „, A 2 ), 

-1 

for every ip £ C X5 (K™ +1 ). Therefore 

(4.64) d T W h + div((l + {k- l)x+(f))VfW h ) - + Aie n , r - A 2 ) 
and, by the definition of Wh, 

(4.65) Wji(- 3 r) = for every r > A 2 . 

Finally by the uniqueness for Cauchy problem [Arj . by (|4.64p and (|4.65|) we 
obtain (gjgg]) . 

Now, performing the change of variable x = h^, t = h 2 r in the integral at 
the left-hand side of j|3~i$jl . we get, by flUEEJ) and (TO!?]) . 

(4.66) /W=ft- n /W. 
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Now, recall that 



(4.67) 



Tc,> [T (;x n ,t; -A 3 e„,0)) 



/4vrf 



In the case k > 1, by Parseval formula, (|3.18j) and (14.67|) we have 



r(i) - 



- m a 2 ,a 3 (C', £; p)F((', -\\\ p)dC,' dx n dtdp 



where Y = M™ _1 x (0, +oo) x (0, A 2 ) x (0, 1) and 
M A2)A3 (C,x n ,t;p) = '\ _ t e — 

(27T)" 1 ^47^ 



x (JC'I + y — k 1 (X " ± A3) VP-) E(C>, x n , A 2 - t; p). 
For fixed A 2 > 0, A3 > 0, taking into account (|3.16p and (|3.17p we have 

M\ 2 .\ 3 (C, x n ,t; p)lm(Ai(p))d(' dx n dtdp 



lim 

Ai^0+ 



Thus 



lim > 



and by (|4.66p the thesis follows. 

Concerning the case < fc < 1, we only give a sketch of the proof, indeed 
such a case can be treated similarly to the case k > 1. In the case < fc < 1 , in 
order to have a suitable formula for first we evaluate the integrals 



Im 



+ OG 



e~ l ^ Xn K(x n + X 3 ,t)dx n 
Im ^ + e-*"*- ^-K(x n + A 3 , 



where 



1 - fc 



v/1 - p|C'| and X(x„,i) 



c -it 



In order to carry out such an evaluation we may use formula 3.322 of |Grj . Then 
we choose A 2 = A?, we perform the change of variable £' = j-, t = Xfr) in the 
integral 1^ and we get 



lim I A? lim ) > 

Ai^0+ \ A 3 ->0+ / 

and the thesis follows. 

Proof of Proposition 13.91 By the triangle inequality we get 

(4.68) |«(i/i,ti;»,t)| = \S 1 (y 1 ,t 1 -y-t)-S 2 (y 1 ,t 1 ;y,t)\ 

> \Sx(yi,tx;y,T)\ - \S 2 (yi,t 1 ;y,t)\. 



□ 
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Let us first estimate from below | £1(2/1, ij.; y, t)\. Recall that T*,(x,t;y,i) is the 
fundamental solution of the adjoint operator of £ + = <9 t — div((l + (fc — l)x+)V). 
Denote by 

Q P /2 = b p/2 (o) x (ti,?), q+ /2 = b+ 2 (o) x (t x ,t), q; /2 = £r /2 (o) x (ti.i). 

By the triangle inequality we have 

(4.69) \S 1 {y 1 ,h;y,T)\>I 1 -R 1 -R 2 , 

where 



(4.70) /1 = 



V X H (x,t;y,t) ■ V x r (x,t;yi,ti)dxdt 



(4.71) i?i = / 



|V x r 2 (x, t; 2/1, t x )| |V x rj (x, t; y, F)|dxdi 



(ti,t)\Q p/2 



Q+ /2 \(D 1 (t 1 ,t)nQ+ /2 ) 



I V x n (x, f; y, f)| | V x r (x, t; y l7 ti)\dxdt, 



(4.72) i? 2 = / |V x r;(x,t;y,t)- V x rt(x,i;y,?)| 

'Di(ti,t)nQ p/a 

x I Vro(ic, j/i, ti)|tixdt 

\V x T (x,t;y 1 ,ti) - V x T 2 (x,t;yi,t 1 )\ 

D 1 (t 1 ,t)nQ p/2 

x\S/ x Tl(x,t;y,t)\dxdt, 

where Di(*i,t) = U te(tl j)fli(t) x {t}. 

Now we estimate from below the term I\ . First we notice that if < 5 < ^= 
then 

(4-73) |x - y| 2 > ^(|x| 2 + p 2 ), \x - yi| 2 > ^(|x| 2 + p 2 ), 
for every x S K." \ -B + /2 - Also, we have trivially 



(4.74) /1 > 



: x(ti,t) 



V x r+(x,i;y,i) ■ V x r (x,t;yi,ti)dxdf 



:+x(ti,t))\Q„ 



|V x rl(x, i; y, i)| |V x r (x, t; yi,ti)\dxdt. 



We now use Proposition 13. 61 and (|4.73[) to estimate from above the second inte- 
gral of the right-hand side of (|4.74[) . We have 



(4-75) 



I^x(t 1> t))\Q+ 



|V x r* (x, t; y, t)\ I V x r (x, t; yi,t\)\dxdt 



p/2 



<c 



s, e 



e c (t-t 1 ) c ft-*) 



" C'o(t-tl) C (t~t) 



;x(ti,t) 



(? - 1)^ (t - h)^ 



R, 
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where Co, Co > 1, depends on k only. Now performing the change of variables 

1/2 



t-tl 



we have 



2p^ 



R<Cq / -= ; /Q ; rfzcfr 



o y/\(l - A) y 

where Cs = max {s n / 2 e _2C, ° s }. By the inequality obtained above, by (|4.74|) 

sg(0,+o o) 

and by Proposition 13.71 we have 

(4-76) h > J— - — , 

where C, C > 1, depends on k only. 

In order to complete the proof we have to estimate from above the terms R\ 
and i?2 defined in (|4.7ip . (|4.72p . Denote by i?n and R\2 the first and the second 
integral at the right-hand side of (I4.71|) respectively. Rn can be estimate in the 
same way of the integral at the left-hand side of (14. 75[) and we have 

(4-77) R 11 < 

where C depends on k only. Concerning R12, by (|2.7a[) and Proposition 13 . 61 we 
have 

ft r ,t) e - Ci (t- H ) e c 4 (t- t) 
(4.78) Rn < C 4 / dt dx' / — p-- -^+rdx n , 

Jt! JR"- 1 J-il>{x>,t) (t - ti)~ (t - t)~ 

whcrc 

iP(x',t) = ~(\xf + \t-t\) 

and C4, C4 > 1, depends on k only. Now we perform, in the integral at the 
right-hand side of (|4.78|) the following change of variables 

t = t 1 +r(t-t 1 ), x' = (r(l - r)fl 2 z\ x n ^{r{\-r)) 1 l' 1 i-\ 1 h. 

Thus, denoting by 

ct(t) = / n v 



nun - II l£ ^===dtdz'dT, 
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we get 

(4.79) R 12 < ^6(h), 
where C depends on k only. Now observing that 

A(z',e,r)>|z'i 2 + (e + ^^vi^) 2 

and applying the Holder inequality, we obtain 

(4.80) 9{h)<C p (±) \ 

for every p G (l,+oo), where C p depends on p and E only. By Q4.77p . (|4.80|) 
and recalling that R\ — Rn + R12 we obtain 

for every h, < h < 5 min{p, v 7 !} and every <5, < (5 < min j ^ , ^ , j , 

where C\ and C2 are defined in (13.23|) and ()3.27j) respectively and C depends 
on k only. 

In order to estimate R2, denote by R21 and R22 the first and the second 
integral at the right-hand side of (|4.72p respectively. By Theorem 14.31 we have 
that there exists a constant C5, C5 > 1, depending on E only such that if 
< 8 < ± and {x,t) G K p := {(x,t) G B p/ c,x{h,t) : x n > c^{\x , \ 2 +\t-t\)} 
then 



(4.82) |V x rj(a:,t;y,t)- V e r;(ar,t;y,t)| < — ■ 



p% (I-t)f+|-f 



where C, C > 1, depends on E and fc only and a — a^j, /3 being defined 
Theorem 14. II We have 

(4.83) R 2 i = J' + J", 
where 

J'= J \V x Tl(x,t;y,T)-V x Tl(x,t;%t)\\V x To(x,t; yi ,h)\dxdt, 
D 1 (t 1 ,t)nQ p/2 nK p 

J" = J \V x T* + (x,t;y,t) ~V x Tl(x,t;y,t)\\y x T (x,t;yi,ti)\dxdt. 

Di(tiJ)nQ p/2 \K p 

By (|4.82p and Lemma [3781 we have 

ha-n 

(4.84) J' < C— , 
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where C depends on E and k only. To estimate from above J" we can arrange the 
method used to estimate i?n and R\ 2 and we obtain that there exists Cg > C5, 
Cq depending on E and k only such that if < 6 < ^ then, for p € (1, +00), 

where C depends on E and fc only and C p depends on p and E only. By choosing 
P = j3a' l F83|) - (USUI and ([4785]) yield 

/-< h a ~ n 

(4-86) fl 21 < - + C— , 

P™ Po 

where C depends on and fc only and C depends on p and i? only. 
Now we estimate R 22 ■ Denote by 

w(x,t) := T (x,t;yi,ti) - T 2 {x,t;yi,t 1 ) 

and recall ()3.28|) . We have that w solves the heat equation in B p / 2 (yi) X 

and, since w(x,t) = for (x,t) £ K™ x (— 00, ti] we can say that w; solves the 

heat equation in B p / 2 (yx) x (t ~ p 2 ,t). On dB p / 2 (yi) x (t — p 2 ,t) we have 

C p 2 C" 

K=M)| < ( t _ tl wa e ~" tt "' 1J X[«i.+°°) - ^Xiu.+oo), 

where C, C" depend on /c only. Therefore by maximum principle and by standard 
regularity estimates we get 

(4.87) |V x w(M)l < in B P/4(yi) x (t~ X>*) • 

It is possible to have a similar estimate for w in (R n \ £? p / 4 (i/i)) x (t — p 2 , t), 

namely by Proposition 13. 61 we have 

(4.88) 

(7 e _ cu*-*i) 

|V x ™(z,i)l < _ < — , (*,*) e (R n \B p/4 ( yi )) x (i-p 2 ,i). 

(i-ii) 2 P 

By (j4~87]) and (14881) we have 

P Jti JR" (t — t)2 p 

where C, C depend on k only. 

The estimate from above of | S 2 (yi ,ti;y,t)\ can be carried out in a similar 
way of that used to estimate the integral in formula (|4.75p . Thus taking into 
account (13.28|) . we get there exists Cj > Cg such that if < S < then 

C 

\S 2 {yi,h;y,t)\ < — 
where C depends on k only. This inequality and (gj35|), fl~56]) . (|4~FT|) . IjiTTS]) . 

dngi, gig give ng^U) . □ 
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